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ABSTRACT 


Random  number  generators  are  notoriously  wasteful  of  digits; 


however,  applyinq  an  augmented  precision  technique  to  a linear  congruential 
generator  enables  one  to  produce  on  even  a small  calculator  a set  of  pseudo- 
random numbers  which  contains  a useful  number  of  elements.  This  paper  sets 
forth  such  a method. 


1.  INTRODUCTION.  Most  modern  computers  and  many  programmable  calcu- 
lators include  in  their  softwares  a function  for  generating  "random"  numbers. 
Such  numbers  are  required  any  time  a "Monte  Carlo"  test  technique  is 
employed. 


It  is  usual  to  tailor  each  algorithm  to  a specific  type  of  use,  and  to  a 
specific  size  of  computer.  Probably  it  is  not  feasible  to  transfer  such  a 
tailored  alqorithm  to  a calculator  of  smaller  si ze--particul arly  to  one  of 
shorter  word  length. 

Perhaps  the  most  efficient  and  certainly  the  most  popular  of  these 
algorithms  is  the  "Linear  Congruential  Generator."  Mathematically  stated, 

x^i  5 (ax^  ♦ c)  mod  m. 

All  quantities  are  considered  to  be  integers.  If  the  modulus  m be  taken  as 
some  power  of  ten  (or  of  two  if  in  binary),  the  modular  operation  is  effected 
by  simple  truncation. 

Most  calculators  have  the  ability  to  truncate  at  the  decimal  point.  A 
decimal  point,  therefore,  is  inserted  solely  for  this  purpose.  Conceptually, 
the  numbers  remain  integers. 
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Given  that  the  modulus  m Is  some  (positive  inteqer)  power  of  ten,  it  is 
found  that  the  algorithm  generates  a full  set  of  m integers  (ranging  from 
zero  to  m-1,  inclusive)  whenever  both 

a = 1 (mod  20)  and 
c -•  ( 1,  3,  7,  or  9)  (mod  lU) . 

The  selection  of  values  for  a and  c is  an  important  part  of  adaptinq  the 
algorithm  to  a soecific  case. 

2.  PSEUDORANOOM  NUMBERS.  Let  us  suppose  that  we  have  defined  a set  of 
m inteqers,  all  different.  A random  selection  from  the  elements  of  this  set 
requires  that  for  any  element,  the  probability  of  selection  be  1/m.  Since 
this  probability  remains  unchanqed  for  subsequent  selections,  sampling  with 
replacement  is  Indicated. 

We  wish  to  develop  an  algorithm  that  does  not  depend  upon  an  outside 
stimulus.  However,  It  remains  necessary  to  provide  a value  for  x0,  so  that 
the  process  can  begin.  This  value  should  be  an  element  of  the  set,  but  the 
choice  can  be  arbitrary.  It  is  called  the  "seed."  After  each  xj  is 
computed  and  used,  it  serves  in  turn  as  the  "seed"  for  the  next  calculation. 
To  avoid  repetition,  some  programmers  employ  a date-time  group  from  which  to 
extract  a value  for  x 0. 

If  any  computed  value  of 


xi  + s ♦ 1 £ (axi  + s + c ) (mod  m) 


7 is  ever  equal  to  some  previously  used  value  of  x-j,  the  algorithm  will  repeat 

itself  over  a subset  of  size  (s+1),  exactly  duplicating  the  previous  cycle. 
If  xj+i  3 x^,  it  is  found  that  s 3 0,  and  the  algorithm  has  already 
degenerated  into  uselessness.  To  circumvent  this,  sampling  without  replace- 


ment Is  used.  But  this  causes  the  probability  of  selection  to  increase  as 


Thus,  the  list  remaininq  element  of  the  subset  can  be  predicted  with 
certainty,  it  is,  of  course,  equal  to  x^,  the  seed  which  began  the  cycle. 

How  then  can  we  presume  to  use  these  sequences  of  numbers  as  "random" 
sequences?  It  is  found  that  if  the  cycle  length  is  very  larqe  (two  hundred- 
fold would  not  be  excessive)  when  compared  with  the  quantity  (of  numbers) 
required,  the  sequence  selected  will  exhibit  certain  of  the  characterist ics 
associated  with  random  sequences. 

The  term  "pseudorandom"  is  used  to  indicate  that  the  sequence  is 
qenerated  by  an  algorithm  so  that  each  element  is  a function  of  its 
predecessor. 

3.  PARAMETER  SELECTION.  At  this  point,  let  us  limit  the  discussion  to 
the  case 

m * 10“ , 

MeM  being  a small,  positive  Inteqer.  Immediately 

v 'irT  * 10®. 

It  was  observed  in  Section  1 that,  under  these  conditions,  maximum  cycle 
lenqth  is  achieved  if  c and  m are  relatively  prime,  and  additionally  a = 1 
(mod  20). 

There  are  other  requirements,  however.  Foremost  among  these  is  the 
restriction  that  ax^  must  never  overflow  the  computer  word  length. 
Should  this  occur,  digits  will  be  lost  from  the  right,  interrupting  the  flow 
of  the  alqorithm  and  seriously  shortening  the  cycle  length. 

The  formula  for  serial  correlation  is 


where  |e|<  4 • 

1 m 

It  can  be  seen  that  the  numerator  varies  from  -0.5  to  + 1.0  and  that  the  two 
terms  are  of  the  same  order  of  magnitude  when 


The  numerator  can  be  reduced  to  zero  by  solving  the  associated  quadratic  in 
c/m.  It  is  found  that 


£ . iii  vo: 

m 2 6 

Now  — v 3 * 0.28867  51345  94812  88225  45/ 43  90  ...  is  irrational,  so  that 

’ 6 

no  element  of  the  set  can  furnish  a value  for  c which  will  reduce  the 
numerator  exactly  to  zero.  It  can,  however,  be  made  quite  small,  whence  "a" 
can  be  set  to  a value  somewhat  less  than  Vm”  without  adversely  affectinq  the 
serial  correlation. 

At  this  point,  it  will  be  instructive  to  examine  the  sequence  generated 
bv  the  following  parameters:* 

x0  * 0 
a * 81 
c * 788677 
m * 1000000 

This  sequence  is  found  in  Table  1-1.  The  entries  are  to  be  read  as  integers. 

It  is  easy  to  observe  that  the  least  significant  diqit  (units  digit)  is 
not  "random"  at  all,  since  it  can  be  predicted  exactly.  In  the  case  at  hand, 


*A11  examples  in  this  paper  will  assume  an  8-diqit  calculator. 


206 


TABLE  1-1 


O . 788677 

0.671514 
0. 181 31 1 
0. 474868 
0.252885 
0. 280462 
0.506099 
0. 782686 
0. 137053 
8 . 939970 

0 . 9 2 6 .?  4 7 
G. 3 14684 
0 . 7 78 0 8 1 
0 . 81 3 2 3 8 
0. 660955 
0.326032 
G. 197269 
0.767466 
0.953423 
0.015940 

0. 079817 
0. 253854 
O. 350851 
0. 207608 
0. 604925 
0. 787602 
0. 584439 
0. 128236 
0. 175793 
O. O27910 

0 - 0 • 1 9 3 3 7 
0. 789024 
0 . 6 9 9 6 2 1 
0 . 457978 
0 . 384395 
0 . 465172 
0 . 467609 
0 . 6 6 5 0 0 6 
0.654163 
0 . 7 5 8 8 0 


CYCLES  OF  DIGITS 


= 81 

» 1,000,000 


0 . 634957 

0. 220194 
0 . 624391 
0 . 364348 
0. 300365 
0. 1 58742 
0 . 646779 
G. 177776 
0. 188533 
0 . 059350 

0 . 636527 
0. 347364 
O . 925 1 6 1 
0 . 7267 1 S 
0 . 652835 
0. 668312 
0- 921949 
0. 466546 
0. 578903 
0. 679820 

0. 854097 
0. 970534 
0. 401931 
0. 345038 
0. 740805 
0. 793382 
0. 0931 19 
0.331316 
0. 625273 
0. 435790 

f. 087667 

0.3  6 9704 
3 . 35470 1 
0. 019453 
© . 364775 
0 . 335452 
0. 960289 
0 . 572086 
0. 1 27643 
0. 127760 


c - h.\,LH 

*o  r 0 


0 . 1 3 7 2 3 7 
0 . 904374 
0 . 083471 
0 . 549828 
0 . 324745 
0 . 093022 
0. 323459 
0. 988856 
0. 886013 
0 . 555730 

0.30  2 8 Q ? 
Q . 816044 
D . 88824 1 
0. 736198 
0. 420715 
0. 366592 
0 . 982629 
0. 381626 
0 . 700383 
0. 5197O0 

0. 884377 
0. 423214 
0 . 0690 1 1 
0. 378563 
0. 452685 
0. 456162 
0. 737799 
0- 550396 
0. 370753 
0.  SI  9670 

0 13  1.  9 4 7 

O. 526384 
0. 425781 
O . 276938 
0 . 220655 
0. 661732 
O.  388  96  _■> 
0. 295 1 66 
G . 697123 
0. 255640 
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it  cycles  through  all  ten  diqits,  then  repeats  itself  exactly.  The  two  least 
significant  digits,  viewed  as  a single  number,  exhibit  a similar  cycle.  A 
good  generator  will  continue  this  effect  until  a cycle  of  length  in  is 
achieved.  If  m is  a power  of  10,  this  maximum  cycle  length  is  obtained 
whenever  both  of  the  following  conditions  are  met: 

1.  a = 1 (mod  20) 

2.  c and  m are  relatively  prime.  This  requires  only 
that  the  final  (ursili)  diqit  of  c be  1,  3,  7,  or  9. 

As  an  aid  to  continuing  the  study  of  the  cycling  effect,  let  us  define 

a$  3 as  (mod  m) 


and 


cs  « c (1  ♦ a ♦ a2  ♦ ...  ♦ as"M(mod  in). 

Given  a * 81,  s * 10,  c * 788677  we  find 

al0  ■ 928801 
cl0  3 939970 

Note  that,  since  x0  * 0,  cl0  appears  in  the  tenth  position  in  Table  1-1. 
Now  cl0  may  be  viewed  as  having  only  five  digits.  It  is  therefore  completely 
exercized  by  a five-digit  multiplier,  and  we  need  merely  use  the  last  five 
digits  of  a10.  The  parameters 


x0  = 0 

3jq  3 28801 

cl0  3 93997 
m j q 3 1 00  000 

will  generate  the  seguence  x l0,  x20,  x 30,  .... 
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At  firct  glance,  this  appears  to  exceed  the  calculator  word  length. 
However,  if  we  multiply  xj  (a1Q-l)  and  then  truncate,  the  algorithm  will  run 
without  difficulty.  To  express  the  complete  formula 


x = x (a  l)(mod  m)  + x ♦ c (mod  m) 

i*10  i 10  i 10 

It  is  convenient  to  compute  as  by  means  of  the  binomial  expansion.  Hence 

( 1 ♦ 80) AO  = 1 ♦ 10(80)  * 45(6400)  + 

♦ 120(80) 3 ♦ 210(80)“  ♦ 

♦ 252(80) 5 ♦ immaterial  terms 

The  previous  strategem  will  thus  be  available  whenever  s is  a multiple  of 
ten.  The  sequence  thus  generated  is  found  in  Table  1-10. 

In  a similar  manner,  the  procedure  can  be  reiterated  and  the  sequence 
xl0Q,  XjQo » x3oo»  * • • generated.  Required  values  of  the  parameters  are: 

*o  * 0 
al00  - 8001 

ci00  - 5197 

ml0Q  * 10  000. 

This  sequence  Is  Illustrated  in  Table  1-100. 

The  process  can  be  carried  no  farther.  To  do  so  results  In  a10oo  * 1*  and 
the  algorithm  degenerates  to  the  successive  multiples  of  xi000.  This  can  be 

observed  by  looking  at  every  tenth  entry  in  Table  1-100.  The  phenomenon  can 

be  called  a "quasi -cycle"  of  length  1000  and  additive  constant  197.  U 
appears  that  original  values  of  "a"  congruent  to  1 (mod  100)  will  hasten  this 
effect  and  therefore  should  be  avoided.  Further  scrutiny  reveals  that  the 
"quasi-cycle"  is  actually  of  length  500  and  additive  constant  598.5. 
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TABLE  1-10. 


CYCLES  OF  DIGITS 


al0  * 28801 

cl0  = MV 

‘ * • • * 
r * - • — • — 

• t 

m = 100,000 

* q =0 

•*.  vv; 

0.93^970 

0. 05377O 

O.i  . i j i 0 

.V’V-.^V* 
, • , * , 1 . 

0 . 01 5940 

0. 574740 

(I.  7 7 7 5 4 0 

•V*  V*V-V“ 

0. 027910 

0 . 02671 0 

O . 6255 1 0 

A A 

0. 775380 

0. 2 14680 

0. 253486 

9 m 9 

0. 059850 

0. 938650 

0.41 7450 

S ‘v*’ 

0. 679320 

0 . 99r*  ;20 

O.917420 

0. 435790 

0. 194590 

0. 553390 

0.  127760 

0. 326560 

0. 125360 

0. 555730 

0. 194530 

0 . 433330 

. » .*  •/ 

0.51 9700 

0.598500 

0 . 277300 

• • 

0.31 9670 

0 . 338470 

0.457270 

;.v.s 

V-jW 

0. 255640 

0 . 2 14440 

0. 773240 

' . ■ «*•  ■ • . 

0. 627610 

0. 026410 

0. 0252 1 0 

. v\ 

0. 735580 

0. 574380 

0.013180 

WA'Xi! 

0. 379550 

0 . 658350 

0. 537150 

• • 

0. 359520 

0. 078320 

0. 397120 

v!v;..v. 

0. 475490 

0. 634290 

0. 393090 

# . * 

0. 527460 

0. 126260 

0. 325060 

• * 

0. 315430 

0. 354230 

0. 993030 

AVy.v 

• • 

0. 639400 

0.11 8200 

0. 197000 

0. 299370 

0.218170 

0. 736970 

0. 095340 

0. 454140 

0.412940 

• • v*/..- 

0. 827310 

0. 6261 10 

0 . 02491 0 

.-  v v -■  #* 
- ■ . * « 

0. 295280 

O . 534080 

0 . 372880 

O. 299250 

0.970050 

f:i . 25605  0 

* .7  • I*  i 

0. 639220 

0. 758020 

0. 476320 

• 

0. 1 15190 

0. 673990 

0 . 832790 

0. 527160 

0. 525960 

0. 124760 

rm  * . i . 

0. 675 1 30 

0.11 3930 

0. 152730 

0 . 359 1 00 

0. 237900 

0.71 6700 

i11  » > 4 

0. 379070 

0 . 697870 

0.61 6670 

•U  --St 

0.  jo'.  04  3 

0 . 293040 

' . 1.:  ' .2(  4 0 

3.62  7'  3 1 0 

G. 8258 1 0 

0. 6246 J 0 

• 

0. 4 5 I960 

0 . 093730 

O. 332380 

V \ 

0.31 8950 

0. 897750 

0. 576550 

■ « ■ » ’ « 

^vv'-V 

0.51 892Q 

0. 037720 

0. 156520 

0. 354890 

0. 313690 

0 . 872490 

• V 

0. 126360 

0. 525660 

0 . 574460 

0 . 634330 

0. 4 73  63  0 

G. 912430 

0. 6 7 S80G 

0. 957600 

0 . 836400 

V ' - • 

* m " ft  * . • . * 

V vv  v ' 

•A-  *\ 

• « 


•V*  -\-V 
> s.- 
*.*/.*■  *•  ,-s 


! !a  .V.'v  „ ';v!v.V  ,lvlvlv5.\ 


TABLE  1-100 


0. 51970O 
G. 639400 
0.359100 
0.678300 
0.598500 
0. 1 18200 
0. 237900 
0.957600 
0.277300 
0. 197000 

0.716700 
0. 336400 
0.556100 
0. 875800 
0.795500 
0. 315200 
0. 434900 
0. 154600 
0. 474300 
0. 394000 

0.913700 
0. 033400 
0. 753100 
0.072800 
0.992500 
0.512200 
0. 631900 
0. 351600 
0. 671300 
0. 591000 

0. 1 10700 
0. 230400 
0.950100 
0. 269800 
0. 189500 
0. 709200 
0:828900 
0.548600 
0. 868300 
0. 788000 


CYCLES  OF  DIGITS 


a 100  = 8001  c i o o = 5197 

mlOO  * 10»000  *o  = 0 


0. 307700 
0. 427400 
0. 147100 
0.466800 
0. 386500 
0.906200 
0.025900 
0. 745600 
0. 065300 
0.985000 

0.504700 
0. 624400 
0. 344100 
0. 663800 
0. 583500 
0. 103200 
0. 222900 
0. 942600 
0.262300 
0. 182000 

0. 701700 
0.821400 
0.541100 
0.860800 
0.780500 
0.300200 
0.419900 
0. 139600 
0.459300 
0. 379000 

0. 8987G0 
0.0184O0 
0. 738100 
0.057800 
0. 977500 
0.497200 
0.616900 
0.336600 
0.656300 
0.576000 


0. 095700 
0.21 5400 
0. 935100 
0. 254800 
0. 1 74500 
0. 694200 
0. 813900 
0. 533600 
0. 853300 
0.773000 

0. 292700 
0. 412400 
0. 132100 
0. 451800 
0. 371500 
0. 891200 
0. 010900 
0. 730600 
0.050300 
0.970000 

0. 489700 
0.609400 
0.329100 
0. 648800 
0. 568500 
0. 088200 
0. 207900 
0. 927600 
0. 247300 
0. 167000 

0. 686700 
0. 806400 
0. 526100 
0. 845800 
0. 765500 
0.285200 
0.404900 
0. 124600 
0.444300 
0. 364000 
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The  conclusion  to  be  drawn  is  this:  Even  though  the  values  of  a and  c 
be  chosen  so  that  the  algorithm  generates  the  full  cycle  of  m integers 
before  repeating,  the  number  of  elements  of  any  "useful"  subset  probably  does 
not  exceed  y ^"m.  What  is  needed  is  a device  to  increase  the  effective  word 
length  of  the  calculator.  How  this  can  be  done  forms  the  subject  matter  of 
the  next  section. 

In  summary,  let  us  view  the  number  axi  + c before  truncation. 
Obviously,  the  left-hand  (most  signifirant)  digits  are  lost  via  the  modular 
operation,  leaving 


xi = ( ax i ♦ c)(mod  m). 

Now  the  xi  can  assume,  at  most,  "m"  different  values.  Therefore,  since  both 
"a"  and  "c"  are  fixed,  the  quantity  axj  * c also  can  assume,  at  most,  "m" 
different  values.  What  this  means  Is  that,  provided  the  values  of  "a"  and  "c" 
are  selected  to  produce  maximum  cycle  length,  the  act  of  truncation  does  not 
reduce  the  quantity  of  numbers--only  their  size.  It  also  shuffles  their 
order . 

What  remains  is,  of  course,  xi+j.  It  is  usual  to  regard  several  of 
the  right-hand  (least  significant)  digits  as  "not  significantly  random."  They 
are  retained,  however,  for  smooth  operation  of  the  algorithm,  and  to  ensure 
that  the  full  complement  of  "m"  different  numbers  is  delivered. 

4.  AUGMENTED  PRECISION  ARITHMETIC. 


Double  precision  arithmetic  is  available  in  the  software  of  many 
computers,  and  even  In  some  calculators.  It  Is  cumbersome  to  program  and 
executes  very  slowly.  This  Is  particularly  true  with  division. 

However,  the  algorithm  for  the  linear  congruential  generator  does  not 
employ  division.  Moreover,  since  a*  < m,  the  word  length  (m  ^m  - 1)  is 
sufficient. 
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Let  m = 102e,  where  e is  any  snail  positive  integer.  The 

"augmented"  word  consists  of  three  parts,  each  of  which  consists  of  "e" 
diqits. 

Let  us  express  x^  in  the  form 

x . ~ u • * 10e  + v • 
i i 1 

Thus  a,  u^,  and  are  all  inteqers  less  than  v/nT,  and  the  product  of  any 
two  of  them  will  not  cause  overflow. 

* 

Some  calculators  will  compute  (but  not  necessarily  display)  an  extra 
riiqit.  For  them,  the  procedure  is  extremely  easy.  First,  compute 

(aui  * 10e)(mod  m) . 

To  this  quantity,  add  (av,  + c)  and  truncate  again.  The  result  is 
Xi+1‘ 


When  place  for  an  extra  diqit  is  lackinq,  it  is  necessary  to  devise  a 
procedure  which  avoids  overflow.  The  following  method,  which  assembles 
x.+1  by  parts,  beqinninq  at  the  right,  works  quite  well. 

As  before,  express  x.  in  the  form 

p 

x.  = u.  * 10  + v j . 

In  analoqous  fashion,  express  "c"  as 

c * p * 10e  + q. 

Store  p,  q,  u^  and  separately.  Select  "a"  so  that 

a < 10e 

a Hi  (mod  20) 
a % 1 (mod  100) 
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18.  Consequently,  multiplication  by 


j ■>. t'.r  f Jit* m n r rr?*  » ' i n«^T *.* i*T 


It  will  be  found  that  a < 10e 
parts  will  net  produce  overflow. 


We  have  immedi ately 

vi  + l H (avi  + Q)  mod  1 0° . 


Since  we  wish  to  retain  both  parts  of  (av  + q),  we  compute  (av  + 


q)  * 10'e,  then  "FRC  ( (av , + q)  * 10'u).* 


-e- 


Vi+i  is  now  stored,  replacing  vj 


Then  u | * [ £ ( a u j + P 


106(avi  + q - v.  + 1))  mod  10e. 

The  sequence  of  nunbers  generated  by 


Xn  = 


a 

c 

m 


0 

9941 

2113  2487 
100  000  000 


is  displayed  in  Table  2-1. 

5.  RANDOM  SELECTION.  RANDOM  ORDERING. 


So  far,  an  algorithm  has  been  developed  which  will  generate  a full  set  of 
m pseudorandom  numbers.  However,  the  length  of  a useful  sequence  of  these 


numbers  is,  at  best,  uncertain  and  doubtless  does  not  exceed  y\/m.’ 


If  a subset  of  far  smaller  but  exactly  known  size  is  to  be  placed  in 
randan  order,  or  if  random  selections  from  its  elements  are  to  be  made,  the 
following  can  be  done. 


*"FRC"  means  "fractional  part  of." 


Storage  space  must  be  provided  to  accommodate  all  the  elements  of  the 
subset,  plus  one  more.  It  may  be  possible  that  scratch-pad  storage  is 

adequate. 

Let  us  illustrate  the  method  by  example.  Suppose  the  task  at  hand  is  to 
shuffle  a pack  of  52  playing  cards,  i.e.,  to  place  them  in  random  order.  We 
thus  require  53  storage  registers,  which  we  number  from  00  to  52,  inclusive. 
The  individual  card  names  are  entered  into  registers  00  through  51  in  any 
arbitrary  order.  N * 52  Is  the  subset  size. 

He  employ  the  generated  sequence  of  numbers  given  in  Table  2-1.  These 
nunbers  (integers)  should  be  distributed  uniformly  on  the  interval  0 to  m. 
Dividing  by  m,  then  multiplying  by  52,  yields  a sequence  uniformly  distri- 
buted on  the  interval  0 to  51.99999  ....  The  "integer"  portion  of  thi: 
nunber  is  used  as  an  address  for  selecting  a card.  That  card  is  then  placed 
in  storage  register  52. 

Next,  all  cards  with  location  nunbers  greater  than  the  "selected" 
location  are  cascaded  downward  one  position.  This  includes  the  card  placed  in 
register  a2.  So  far,  the  illustrative  example  has  given  52  x 0.21132467  = 
10.988  ....  The  card  in  location  10  was  drawn  and  stored  in  location  52. 
Say  it  is  the  Spade  Jack. 

After  cascading,  only  51  cards  are  of  Interest.  Hence  51  * 0.99185754  = 
50.584  ....  The  card  in  location  5Q--the  King  of  Clubs--is  drawn  and 
placed  in  register  52.  Again  after  cascading,  the  subset  of  unshuffled  cards 
is  reduced  to  50  in  nunber.  Hence  50  x 0.26713001  * 13.356  ....  The  card 
now  in  location  13--the  deuce  of  Hearts--is  drawn  and  placed  in  register  52. 

Continuing  as  above,  49  x 0.75075428  * 36.786  ....  Tht  card  in 
location  36--say  the  King  of  Dia»onds--is  selected  and  placed  in  location  52. 

When  the  size  of  the  unshuffled  subset  is  reduced  to  unity,  that  card 
certainly  will  be  found  in  location  U0,  and  it  certainly  will  be  selected  for 
transfer  to  location  52.  Consequently,  that  transfer  can  be  effected  without 
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TABLE  2-1 
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CYCLES  OF  DIGITS 


a «=  9941  c = 2113  2487 

m = 100,000,000  x0  = 0 
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0.211 324*7 
0.99185754 
0. 267130O1 
0. 75075428 
0. 45962235 
0. 31710622 
O. 56425789 
0. 49900936 
0.86337263 
0. 9986397© 

0. 68858257 

0.41O65324 
0.51518371 
0 f 65258598 
0 . 56855205 
0. 18725392 
0.70254359 
0. 19715306 
0. 10989433 
0. 67085940 

0. 22462027 
0. 16142394 
0. 9764174 1 
0. 77679768 
0. 35706175 
0. 762181 62 
0. 05380929 
0. 33447676 
0. 74479603 
O. 228659 1 G 

0.31 143797 
0. 21618464 
0.'  30283  111 
0.65538938 
0.43715145 
0. 93388932 
0. 00505499 
0. 46293046 
0. 70007773 
0. 68403830 


0.241 03567 
0. 34692034 
0. 9 4 6 4 2 4 3 i 
0.6203610  8 
0. 220821 1 5 
0. 39437702 
0 . 71328 u 6 9 
0 . 9 3 4 6 1‘  4 1 6 
0. 70773943 
0 . 84899850 

0 . 10541337 
0 . 12563604 
0. 15919851 
0 . 8037 1 278 
0. 92007O85 
0.63564472 
0. 15548639 
0.90152786 
0.299781 1 3 
0. 33553820 

0. 79657107 
0. 92433174 
0. 9931522 1 
0. 1 3744448 
0 . 54690055 
0. 94969242 
0. 1 0367209 
0.81557156 
0. 80820283 
0 . 55565798 

0. 00650877 
0. 91500744 
0 . 3002859 1 
0. 35355618 
0.91331 025 
0. 428520 1 2 
0.  1 2983779 
0 . 92879526 
0 . 3650O453 
0.721 35760 


U . 22722647 
0. 06966314 
0. 73259961 
0 . 9y4  04  i y y 
0.63129995 
6. 964 12782 
0 . 60598349 
O. 29319896 
0. 90218623 
0. 84463730 

0. 750724 1 7 
0. 16029884 
0. 74209331 
0. 36091958 
0. 11286965 
0.24851552 
0. 70410919 
0. 76073266 
0. 15174793 
0. 73749700 

0. 66900187 
0. 75891454 
0 . 58076701 
0.61617128 
0. 57001935 
0. 77368322 
0. 39621489 
0 . 98354636 
0 . 64568963 
0.01 193670 

0. 87405957 
0. 23751024 
0 . 3006287 1 
0.  68180298 
0. 01474905 
0 . 83 1 63092 
0. 45430959 
0.41 349806 
0.71601133 
G.  0799564G 
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employing  the  algorithm.  Further,  cascading  can  be  omitted  or  not,  at  tne 
pleasure  of  'he  programmer. 


The  final  result  is  the  shuffled  deck,  in  order  of  selection,  in  the 
des.  ;ated  storage  locations.  Omitting  the  final  cascading,  the  example 
leaves  the  Spade  Jack  in  01,  the  Club  King  in  02,  the  Heart  deuce  in  03,  the 
Diamond  King  in  04,  etc.  The  shuffled  deck  can  new  be  put  to  the  use  for 
which  it  was  intended. 


If  there  is  a requirement  to  "deal"  the  cards  one  at  a time,  it  is 
suggested  that  the  card  in  the  highest  numbered  location  be  taken  first.  Not 
only  is  the  prograinmi  ng  simpler,  but  the  stigma  is  avoided  which  usually  is 
attached  to  dealing  from  the  bottom. 


In  summary,  a set  of  uncertain  size  has  been  used  to  produce  a much 
smaller  subset  of  known,  fixed  size. 

6.  STATISTICAL  TESTS.  There  is  much  to  be  found  in  the  literature  on 
the  subject  of  testing  sequences  of  numbers  to  determine  whether  or  not  a 
sequence  could  have  been  produced  by  a random  selection  process.  These 
methods  will  not  be  repeated  here. 


It  is  enough  to  be  reminded  that  the  answers  to  these  statistical  tests 
will  be  stated  as  probabilities.  We  should  read  nothing  into  the  result 
beyond  the  probability  statement  itself. 
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